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Abstract 

Transport properties of one-dimensional Kronig-Penney models with binary cor- 
related disorder are analyzed using an approach based on classical Hamiltonian maps. 
In this method, extended states correspond to bound trajectories in the phase space 
of a parametrically excited linear oscillator, while the on site-potential of the original 
model is transformed to an external force. We show that in this representation the 
two probe conductance takes a simple geometrical form in terms of evolution areas 
in phase-space. We also analyze the case of a general iV-mer model. 
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1 Introduction 



The random dimer model, a tight binding model with correlated disorder, introduced in 
refs. [|], |2| has attracted considerable attention due to the presence of transparent states 
in an otherwise disordered one dimensional system |J-[(|]. In the present paper we address 



the issue of the spectrum of the random dimer Kronig-Penney (RDKP) model Jl0| and 
extensions using the same Hamiltonian approach applied earlier in the context of the tight 
binding models ||. Using techniques from dynamical systems theory [|11|], we construct a 
Poincare map that turns the Kronig-Penney model into an equivalent tight-binding model 
and study the latter through a two-dimensional map corresponding to a classical linear 
oscillator with a parametric perturbation given in the form of periodic delta- kicks [f| . The 
amplitudes of these kicks are defined by the site potential of the tight-binding model. In this 
representation, extended states of the tight-binding model are represented through bounded 
trajectories in the phase space of the Hamiltonian map. Furthermore, in this representation 
the two probe conductance is related to the time evolution (under the Hamiltonian map) 
of areas initially defined by the basis unit vectors. This new approach provides an effective 
and simple tool for understanding transport properties, the structure of eigenstates as 
well as for deriving analytical expressions. In particular, one can easily determine fully 
transparent states for the general case when N -sites are correlated. In the following 
section we summarize briefly the Hamiltonian map approach used in ref. [0] and apply it 
in the context of the RDKP model. In section 3 we analyze transport properties through 
a new expression for conductance while in section 4 we conclude. 

2 Hamiltonian Map Approach 
2.1 Time-Dependent Linear Map 

The model of interest is the one-dimensional Schrodinger equation with an array of 5— function 
potentials: 

E<t>{z) = + £ e n S(z - z n )<t>(z). (1) 

a Z n=l 



Equation (0) defines the Kronig- Penney model where E is the eigenenergy of the stationary 
states, e n is the strength of the potential and z denotes space while we take the positions 
of the S functions to be regularly spaced (z n = n). The tight-binding model corresponding 
to Eq. (@) is |1 



S1TI ( Q ) 

<Pn+l + 0n-l = V n <f) n , V n = 2cOs(q) + 6 n (2) 

Q 

with q 2 = E, (j) n = <ft(z = n); Eq. (||) can be written equivalently as a two-dimensional 
map M n , i.e.: 

( ) - ( ? V )(::)• < 3 » 

where 4> n = x n and y n = n -i- An eigenstate of Eq. (0) is a 'trajectory' of the map 
of Eq. (13]). Straightforward diagonalization of this map leads to the eigenvalues = 



_ e ± lfln ^gj-g ph ase ^ j g introduced by the relation v n = 2 cosfi n . For 
\v n \ < 2 we obtain stable map rotation with phase /i n . For \v n \ = 2 we find the curves 
separating regions of allowed and forbidden energies, determined through equations q = 
(2k + l)vr, for e = -2qcot(q/2) and q = 2kir, for e = 2qtan(q/2) with k = 0, 1, 2, . . .. To 
understand the origin of the resonant states in the RDKP model resulting when adjacent 



pairs of random energies e n coincide [|1(| , we consider the sequence e n which consists of one 
dimer only i.e. where all the values of e n are equal to ex except two values for which we have 
Cm = = £2- From the property of the map eigenvalues found previously we observe 
that this unique dimer with energy €2 does not influence the trajectories of the map of Eq. 
d^D when the total phase advance \i m + fi m+ i = 2/i m through the dimer is equal to n or 
2tt . Since the latter value /i m = 7r is forbidden from the stability conditions, the resonant 
energy q 2 cr = E cr is defined by /i m = %/2 giving —2q cr /e 2 = tan(q cr ). As a result, for the 
general case of randomly distributed dimers e\ and e 2 , there are two resonant values, 

e 12 tan(q cr ) 

Qcr = (4) 

for which dimers of the first type, e\ , or of the second, €2 , have no influence on the 
transparent states. One should note that since tan(q cr ) is a tt— periodic function and takes 
values in (—00, +00), we have an infinite set of critical energies q cr (two in every interval 
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[(2k — l)7r/2, (2k + 1)tt/2}) flQfl . However, for the first allowed band defined through the 
condition \v n \ = 2, the disorder strength e must be greater than a critical value e cr = —2 
in order to have resonant states. The critical disorder for reflectionless modes to appear, 
arises only in the first band of the spectrum. 

The previous analysis can be extended to the general case of iV-mer (two values e± 
and €2 appear in blocks of length N) where the resonant energy is defined through the 
following condition 

7r 2n 3n 4n (j + 1W „ T . , 

We note that as the block size iV increases the number of resonant states proportional to 
2(N — 1) increases as well. In the first zone, in particular, the disorder strength should be 
smaller than a critical strength e cr (N) for this to happen. The latter is obtained through 
equation 2cos(q cr ) + e 1)2 sm ^ q ° r ^ = 2 cos fiN that gives the corresponding critical wavevectors. 
Equivalently, e cr is obtained through: 

sinq 2q 

= • (6) 

cosq — cosun e cr 

The resonances inside the first band appear whenever the derivative of the left hand side 
of Eq.(|) at q = is less than — j- and thus 

e cr (N) = 2(cos(fi N ) - 1) (7) 

In the phase diagram of Fig. 1 we distinguish three regions: (I) contains no resonant states, 
(II) where some resonant states appear and (III) where all resonants are concentrated. The 
borders of region (II) start at e = —2 for iV = —2 (RDKP case) and are given by the curve 
Eq. (0) with fi N = tt/N (upper bound) and //jy = (N — l)n/N (lower bound) and approach 
zero and —4 respectively as iV — ► 00. We recall that for the perfect Kronig- Penney lattice 
the range of accessible e values within the first energy band is [—4, +00). We further note 
that the lower border corresponds to the existence of only one resonant state, while as 
we are going towards the upper border more resonances appear. The upper critical curve 
delimits region III in which all system resonances are found. 



By introducing a new variable p n +i = Xn+i — %n playing a role similar to momentum 
we obtain through Eq. (|]) a new map representation: 

Pn+1 Pn fn -^n ^g^ 
•En+1 Pn+1 

where /„ = v „ — 2 and having the same eigenvalues with the original map (D) . In the map 
of Eq. (H) ellipses correspond to e n = e\ for all n, a defect at site m with e m = e 2 results in 
a kick into another ellipse (Fig. 2a) while correlated defects lead to a return to the original 
ellipse, since the total phase advance /i = fi m + /i m+ i = n (Fig. 2b). When a random 
mixture of dimers with energy e 2 is embedded in a chain with energy ei, we obtain a phase 
space trajectory similar to the one in Fig. 2c. We observe two ellipses corresponding to e± 
and e 2 values respectively. The second ellipse is formed by points occurring every time the 
first site of a dimer is encountered by the map. In the general case of correlated M— block 
under the condition that the total sum J2n=i A*n of phase shifts is equal to mn, for any 
sequence with m integer, the trajectory always returns to the ellipse associated with the 
"perfect" sites to the left and to the right of the scattering potential. 

2.2 Parametric Linear Oscillator 

Another useful representation of the original model of Eq. (|2|) similar to the map of Eq. (^) 
but more convenient for the analysis of the localization length, can be obtained through 
two successive maps ||: 

Pn Pn A n X n , , 

•^n *^"n. 



and 



Pn+1 =Pn COS /X ~ X n Sm/I , , 

x n +i = Pn sin Ho + X n COS/i 



When the maps of Eq. (||) and Eq. ( |T0| ) are combined, they result in a form of Eq. (0), 
viz. 

x n+ i + x n _i = (2cosyU + v4 n sin/i ) x n . (11) 

Comparing with Eq. (H) one can establish the correspondence fio = q; A n = y 
between the parameters q, e n in the original model (0) and the parameters /iq, A n of the 



map of Eq. (|9|, [T0|) . The latter map has a clear meaning since the map of Eq. (g) corresponds 
to an instant linear kick of the strength A n resulting in the change of the momentum p n 
and the map of Eq. (|ToD describes free rotation in the phase plane (p, x) defined by the 
angle Pq. The dynamical system modeled by Eqs. * s ^ na ^ °f a li near oscillator with 

a periodic parametrical perturbation with a Hamiltonian || 

H = + ~\x%{t) ■ Ut) = £ A n 6(t-n) (12) 

n=— oo 

We note that by integrating Eq. (|l|) between two successive 5-kicks of the potential, we 
obtain Eq. (|10D , while integration over a kick, results in Eq. (||) with a new kick strength 
A n — e n . By comparing Eq. (JJ) and Eqs fl9|),(|T0D we find the significance of the variable p n : 
it is the rescaled (with respect to q = \/~E) first derivative of the local amplitude function 
4> just before the n— th kick i.e. p n = (d<p/dz) z=Zn /q. 

For the dimer case defined by the two values of 61,62, we can set without loss of 
generality e\ — 0. As a result, the motion corresponding to e n = e.\ is represented by the 
circle in the phase plane (p, x) and resonant behavior results when, after a given number 
of kicks with e n = 62, the trajectory returns to this circle. An example of this behaviour is 
given in Fig. 3a for q = q cr = 11.4 and €2 = 9.75885. A similar behaviour for the case of a 
trimer = 3 is illustrated in Fig. 3b for 62 = 1.25 and q = q cr = 1.489. We observe that 
the trajectory is bounded in the phase-space. 

The above scheme is also valid for the much more general case where the locations 
of the 5— functions z n in the original model Eq. (|]) are not equidistant and can be taken 
from an arbitrary distribution. In this case Eq. ( |i~0D has the same form provided we make 
the substitution /x — > p>T n — QT n , where T n corresponds to random periods of the kicks in 
the Hamiltonian approach (0). In the original model of Eq. (j^), T n indicates the random 
distance between two successive lattice sites, i.e. T n = z n+ \ — z n . The critical value q cr 
for the dimer is obtained directly from Eq. ( |TTD using the relation 2cos/xat = 2cospx n + 
^4 n sinyUT n with p,^ = n/2 (see Eq. (||)). As a result, we obtain tan(g cr T n ) = — e/(2q cr ) 
where both T n and e n are related to the same lattice site z n in Eq. (0). Therefore, we 
conclude that in the case of a generalized dimer where time displacements T n and on-site 
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potential are paired, in a way such that T n = T n+ i = T 2 and e n = e n+ \ = the condition 
for the critical energy will be similar to the one obtained previously but with the change 
to -> HT n =fMi = qT 2 . 

2.3 Nearly Resonant States 

The representation of the model of Eq. (Q) used in section 2.2 allows for the study of global 
properties of eigenstates. In particular, the resonant delocalized states correspond to a 
bounded motion described by the maps of Eqs (PD,(P),([10|). Localized states on the other 
hand are represented by unbounded trajectories that escape from the origin of phase space 
(p, x). This is illustrated in Fig. (4) for the case of random dimers with non- resonant values 
of q . The exponential increase of a distance from the origin (p = x = 0) is related to the 
localization length of the corresponding eigenstate. In order to study the dependence of the 
localization length I for nearly-resonant states, it is useful to pass to action-angle variables 
(r,9) for the map of Eq. ((|| (|i~0|) using the definitions x = r cos 9 and p = r sm9. We 
obtain a map for the action r given by 

rl+i = r 2 n D 2 n , D 2 n = (l + A 2 n cos 2 9 n + A n sin 29 n ) (13) 
where the transformation for cos 9 n and sin 9 n is given by the relations 



cos 9 n+ i = D n 1 {cos(6' n + /i ) - A n cos 9 n sin /i } 
sin0/ n+ i = D~ l {sin(9 ri + ji ) + A n cos0„ cosyU } 



(14) 



The relations of Eq. (|T3"|,I4") can be used instead of the common transfer matrix approach 
for the determination of the localization length I . The latter is equal to the inverse of the 
Lyapunov exponent 7 defined as 

1 N ~ 1 r , 

7= lim V ln(— ) (15) 

where the ratio r n+1 /r n = D n is given by (|T3|). 

The advantage of this approach in the finding of the Lyapunov exponent 7 in compar- 
ison to the standard transfer matrix method is that there is no divergence during iterations. 
It is interesting to note that Eqs. (|i~3D , (|i~4l) can be mapped into a one-dimensional map 



9 n+ i = F(6 n ) which is non-linear for the non-zero perturbation A n ^ 0. One can show 
that in such a representation, the expression fllTf ) is directly related to the stretching of the 
phase, d6n+i/dO n = D 2 . Therefore, the original quantum problem is reduced to the study 
of the properties of a one-dimensional time-dependent map and its tangent space. 

Due to correlations in the sequence 9 n , the expression (|T5|) can not be evaluated 
directly However, it is possible to construct an effective map for two successive kicks 



of the single map ([13]) and neglect the correlations between the phases 6 n+ 2 and 6 n near 
the resonance q = q cr — 5 ~ q cr . Applying the resulting two-step map to the present 
case allows us to estimate the Lyapunov exponent for 8 -C 1, using the expansion in 
W = A n 5 1 sin/x , with the successive averaging over 6 n , leading to 

i^Q^r- ( 16 ) 

The factor Q stands for the probability for the dimer of the second kind (with energy €2), 
to appear. In Fig. 5 we compare the analytical result (|I~6"D (solid line) with the numerical 
data obtained from the map (|T3|JI^) (circles) after iterating up to 4000000 time steps and 
averaging over more than 1000 realizations of the random kick strengths A n , for the case 
Q = 0.5, e = 9.75885 and q cr = 11.4. The agreement between theory and numerical data is 
extremely good as seen in Fig. 5 verifying that expression ([16]) valid for 5^1. 

From Eq. flTB"]) we determine the dependence of the inverse of the localization length 
for the near-resonant states and the way it changes when the system parameters change. 
We note that the higher the order of the resonance, q » 1 (i.e. the higher k in [(2k — 
l)7r/2, (2k + l)7r/2}), the larger is the localization length, and thus better transport prop- 
erties are expected. Such a behaviour of l(q) is expected since for q » 1 the second 
term of v n (where q appears in the denominator) in Eq. (^]) becomes negligible resulting 
to a tight-binding equation with zero on-site potential. The localization length increases 
also when decreasing e towards zero. This is easily comprehended since when e\ = £2 = 
we recover the properties of the perfect lattice. Finally, as the concentration Q of dimers 
decreases, the value of the localization length for the near resonant states increases. 
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3 Transport Properties. 



In this section we examine the transport properties of our system by studying the behavior 
of the transmission coefficient. We assume that the system of Eq. (0) is a sample consisting 
of L lattice points with two identical semi-infinite perfect leads on either side. As a result, 
the left lead extends from — oo < n < 0, the sample extends from 1 < n < L and the right 
lead extends from L + 1 < n < oo. The purpose of these leads is to carry the incoming, the 
reflected and the transmitted waves. Here E (see Eq.(|l|)) is the Fermi energy and without 
any loss of generality we choose e n = everywhere in the leads. 

In order to calculate the transmission amplitude ti of a segment containing L sites 
we inject a particle from — oo with an energy E = 2cosq towards the sample. While the 
particle passes through the sample it undergoes multiple elastic scattering. Eventually, it 



we 



comes out of the sample from the right end with amplitude t^. Following Pichard [|12 
write the transmition coefficient Tl = \ti\ 2 in terms of the matrix elements of the total 
transfer matrix P L = lln=i M n as 

4\sinq\ 2 

L = \(P L hi - (Pl)i2 + {Pl)^ - (PL)ne-^r ( } 

In the Hamiltonian map approach, the above system corresponds to the parametric 
linear oscillator of section 2.2 where the strength A n of the instant linear kick (see Eq. (|9])) 
is equal to zero for times i<0ort>L + l describing free rotations in phase plane, while 
in the time interval 1 < t < L the strength A n is determined by the disordered site energy 
e n of the underlying one-dimensional Schrodinger equation ([[]). 

In order to establish a relation for the transmition coefficient in the frame of our 
Hamiltonian map approach, we recast the two successive maps of Eqs (p|) , (|10D to the fol- 
lowing two-dimensional map Q n 



x n+ i \ _ ( cosfio + A n sinfi sin^o \ / x 
Pn+i j I A n cosfi - sin/io cos^ )\p 



(18) 

which is related to the transfer matrix M n defined in Eq. (|3p through a similarity transfor- 
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mat ion R 

Q n = RM n R~ l ; R = ( eoL, J (19) 

\ sirifiQ sinfio / 

From the above Eq.(|TT?D and Eg. (|T7|) , one obtains for the transmition coefficient Tl of a 
system with L scatterers: 

4 

Tl = mE + (Fl%) + {{F L % + (Fl) 22) + 2 (20) 

where the matrix Fl is the product transfer matrix i.e Fl = Iln=i Qn- From Eq. (pO|) 



we see that the sum inside the first parenthesis in the denominator is equal to the inner 

1 




product of the vector v(t = L) = Fl ( n ] i.e. to the modulus square of the vector 



v{0) = (J) evolved under the dynamical map fllSD (or equivalently under the map (||, 
[10|) ) for time t = L, in the phase space of the parametric linear oscillator described by the 
Hamiltonian (0). Similarly, the sum inside the second parenthesis in the denominator in 
Eg. (p0|) corresponds to the modulus square after the evolution for time t = L of the initial 
vector u(0) = ^ J V It is interesting to note that the initial vectors v(0),u(0) correspond 
to the unit vectors pointing to the two perpendicular directions on the phase-space plain. 

Using these observations, we can give a geometrical interpretation for Eq. (|20D viz. to 



relate Tl with areas in the phase-space of the two dimensional Hamiltonian map ([12]). In 
particular, we can interpret the sum inside each parenthesis in the denominator, as the area 
of a circle described by a radius rx, 7*2 which is given by the time evolution (under the map 
dl8|) ) of the initial vectors (ri^, #1,2)4=0 = (1, 0), (l,7r/2). Thus Eq.(pO|) can be rewritten 
in the following form 

T ^ 2S l* (21) 

where S^ ot = S*f + S 1 ^ is the sum of the areas defined by the radius rx, r 2 at time t = L. In 
the case of perfect lattice where we have simple rotations of the initial vectors (ri^, #1,2)°) 
the areas defined after time L will be the same S^ ot = S® ot = 2n and hence Tl — 1. 

Using Eq. (|T3|) we can write Eq. (|21"D in a way that it is more tractable for numerical 
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calculations, i.e. 



T L = 5 5 (22) 

where D l=1 ' 2 correspond to initial conditions (r^,^) = (1,0), (l,7r/2) respectively. 

The results we have obtained so far provide an exact, although non-closed, analytical 
description of any one-dimensional system that can be written in the tight-binding form of 
Eq.(Q). We will now evaluate them for the specific case of RDKP to describe those relevant 
features of the transmission coefficient that may be the fingerprint of extended states. 

In Fig. 6a,b we show numerical results for a system of 10000 scatterers after averaging 
over more than 10000 different realizations of the disordered sample. We take the values 
of ei = and e<i = 9.75885 and the defect concentration Q = 0.5, i.e. the most random 
case. In Fig. 6a we used the critical energy q cr = 5.443223 lying inside the second zone 
of the spectrum while Fig. 6b q cr = 33.13 lying inside the 10 th zone. We note that states 
close to the resonant energies have very good transmission properties, similar to those at the 
resonant energy where the transmition coefficient T is equal to one. This is compatible with 
the findings of the random dimer model with one-band (tight-binding approximation) ||. 
Moreover in RDKP model, the width of the peaks depends on the order of the resonance, 
as mentioned previously (see also |TIJ). From the comparison between Fig. 6a and Fig. 6b 
we see that the higher the resonance the wider the band of states with T « 1. In Fig. 6c we 
presented results for a different defect concentration Q in order to study the dependence of 
the transmition coefficient depends on Q. We use the same values of ex, €2 and Q = 0.2. By 
comparison with Fig. 6a we conclude that as Q decreases the number of transparent states, 
i.e states with transmition coeffissient close to one T ~ 1, increases, in perfect agreement 
with the results of the previous section for the localization length of nearly-resonant states. 

4 Conclusions 

We have studied a Kronig-Penney model with binary on-site disorder randomly assigned 
every two sites. For such a model it was found [|T(J that there exist infinite number of 
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special energies E cr = q^ r at which transparent states appear. We recover these results 
using a new approach based on classical Hamiltonian maps. We generalized our results 
for a TV-mer case and obtained a simple expression for the resonant energy values. We 
constructed a phase diagram in the e cr — N plane which marks three distinct regions: one 
where no resonant states appear, the second that contains some of the resonant states and 
the third that contains all possible resonant states for different correlated blocks of size 
N. This separation in three distinct regions is valid only in the first zone and as a result 
it might have some relevance to the low temperature system properties. Our dynamical 
system approach maps resonant delocalized states to bounded trajectories, while localized 
states are represented by unbounded trajectories in the phase space (p,x). Making use of 
an expansion in the vicinity of the resonance, we derived an analytical expression of Eq. 
fllBD for the Lyapunov exponent for the nearly resonant states. Finally, in the frame of our 
Hamiltonian map we established a simple geometrical picture for the transmition coefficient 
showing that it corresponds to evolution of areas in the phase space of a linear parametric 
oscillator. Using these last results we calculated the transmission coefficient that exhibits 
peaks up to T = 1 for energy values equal to the resonant ones. Near the resonant energies 
there are nearly-transparent states with large transmition coefficient, the number of which 
is inversely proportional to the defect concentration Q and increases with the resonance 
order. The properties of the RDKP model that were analyzed in this work could be used 
to several mesoscopic quasi-one dimensional studies. 
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Figure Captions 



Figure 1: The phase diagram showing the critical e cr values as a function of the block size 
N. 

Figure 2: Phase space of the map of Eq. (13) for p = x = 1: 

(a) One value of e 2 in the sequence e :...e\t\t\t2t\t\t\... for q cr = 11.4; 

e i = 0;e 2 = 9.75885 ; We note that one point is outside of the ellipse representing 
the kick to the €\ trajectory by e 2 . 

(b) Two values of e 2 (one dimer) in the sequence e : ...e 1 e 1 e 1 e 2 e2eiei e i--- for e i = ; 
e 2 = 9.75885; q cr = 11.4. We note that there is a point inside the ellipse representing the 
kick to the e\ trajectory by the first e 2 value. The second e 2 value kicks the trajectory back 
to the ellipse. 

(c) Dimers of type e 2 , randomly (with probability Q = 0.5) distributed in the se- 
quence e for ei = 0, e 2 = 9.75885 and q cr = 11.4. 

Figure 3: The phase space of the map (16-17) for p = x = 1 and Q = 0.5, e\ — : 

(a) N = 2 ( dimer); q cr = 11.4; e 2 = 9.77885; 

(b) iV = 3 (trimer); g cr = 1.489; e 2 = 1.25; 
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The length of sequence e is equal to L — 1000 . 

Figure 4: Nearly resonant states for a dimer (N = 2). Comparing with the resonant states 
shown in Figs.2-3a, nearly-resonant states correspond to the unbounded (for t — > oo) mo- 
tion with a slow spread of the points in the phase space: 

(a) q = 11.399 ; e 2 = 9.75885; 

(b) q = 11.390 ; e 2 = 9.75885; 

(c) q = 11.38 ; e 2 = 9.75885; 

Figure 5: Numerical (circles) and analytical estimation (solid line) of the Lyapunov exp. 
for near resonant states. We use q cr = 11.4. We observe good agreement. 

Figure 6: Transmission coefficient (averaged over more than 10000 realizations) as a func- 
tion of q = E 2 for a system with e 1 = 0, e 2 = 9.75885: 

(a) Q = 0.5 and q cr = 5.44 corresponding to a resonance inside the second band. 

(b) Q = 0.5 and q cr = 33.13 corresponding to a resonance inside the 10 th energy band. 
We note that with respect to (a) the band of states with T 1 is wider. 

(c) Q = 0.2 and q cr = 5.44 corresponding to a resonance inside the second energy 
band. We note that with respect to (a) we have now more states with T« 1. 
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